Enhanced mesoscopic correlations in dynamic speckle patterns 
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We study the temporal evolution of speckle patterns in transmission of short wave pulses through 
a disordered waveguide. In the diffuse regime, the short-range spatial structure of speckles is the 
same as for continuous wave (cw) illumination, whereas the long-range correlation between distant 
speckle spots grows linearly with time and can exceed its cw value. We discuss the physical origin of 
this phenomenon, compare our results to recent microwave experiments, and suggest that a similar 
linear growth with time should also be characteristic for other mesoscopic interference phenomena. 
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In a disordered medium of size L, waves can either 
propagate ballistically (if L <C mean free path £), by 
diffusion (if L ^ £), oi can be localized (if L exceeds 
the so-called localization length ^ > £) QlEllli. The 
different regimes of wave propagation can be probed by 
sending into the medium a short wave pulse and by study- 
ing the decay of the intensity I{t) of waves transmitted 
through the sample. In the ballistic regime, I{t) mimics 
the incident pulse, whereas in the diffuse regime I{t) ex- 
hibits appreciable sample-to-sample fluctuations with an 
ensemble average (/(i)) that decays exponentially with 
time l^, |5j- The characteristic time tjj ^ L^/D of this 
decay (the Thouless time) features the diffusion constant 
D. In the localized regime, the decay of (/(i)) slows 
down, suggesting that D ^ |2,l3l- A similar situation 
is encountered when a brief voltage pulse is applied to a 
disordered conductor (with I{t) staying for the current 
through the sample) [ajl3ila| or when a quantum particle 
is injected into an open chaotic cavity (with I{t) staying 
for the survival probability) j2| . 

Recently, it has been realized that already in the diffuse 
regime the closeness of the localization transition causes 
the average intensity (/(i)) to deviate from pure exponen- 
tial decay [3, [l^ |lj ll3| • More precisely, the apparent 
diffusion constant decreases linearly with time t until t 
becomes of the order of the Heisenberg time tn ^ to- 
Beyond the Heisenberg time, the so-called anomalously 
localized states dominate the dynamics p LD, la]- Thus, 
as time goes on, (/(i)) behaves as if the system were un- 
dergoing a continuous transition from the diffuse regime 
for i ~ iu to the localized regime for t > tn, when 
extended states have died out and only localized states 
are relevant. However, the intensity of waves in a disor- 
dered medium fluctuates rapidly in space |12i] , giving rise 
to the so-called speckle pattern, and hence the average 
value is not sufficient for its statistical description. The 
speckle pattern is observable by eye for light {see^e.g., 
Ref. [lj|) and readily measurable for microwaves [la| or 
ultrasound [la|- Does the speckle pattern /(r, i) evolve 
with time and if it does, in which way? Answering this 
question constitutes the principal subject of the present 
work. 

The speckle pattern can be quantitatively character- 



ized by the correlation function of intensity fluctuations, 
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where 51{r,t) = /(r,t) - (/(r,t)). Up to now, C(Ar, i) 
has only been studied for continuous-wave (cw) illumina- 
tion in which case it is independent of t. C can be split 
into a sum of short- and long-range 'mesoscopic' con- 
tributions (ij, [13 ■ The short-range contribution (usu- 
ally denoted by Ci) is of order I for Ar < wavelength 
A but decays rapidly for Ar > A, identifying A as the 
size of an individual speckle spot [ij, [iSj . Much atten- 
tion has been given to the long-range part of C (denoted 
by C2), which survives for Ar > £ and can be inter- 
preted as a correlation between distant speckle spots 
|l5l llSl I19I I2OJ. Detailed theoretical and experimen- 
tal studies of C2 have been undertaken for waves in 
a quasi- ID disordered waveguide jl^ |2flj, and we will 
also limit our consideration to this case. In the dif- 
fuse regime (dimensionless conductance of the waveguide 
g 3> 1) and for cw illumination, it has been shown that 
C2(Ar,i) = C^""'^(Ar) ~ I/5 < 1 [3 El 0- For 
pulsed illumination, the evolution of C(Ar, t) with time 
has been recentl y st udied in microwave experiments by 
Chabanov et al. |2]| . The experiments suggest that Ci is 
independent of t and that C2 increases with time. How- 
ever, no theoretical model has been proposed so far to 
explain these observations. 

In the present Letter we show that when a pulse of du- 
ration tp is sent into a disordered waveguide, the short- 
range spatial structure of speckle pattern measured in 
transmission through the waveguide is essentially the 
same as for cw illumination (speckle spot size ^ A). In 
contrast, the long-range structure of speckles appears to 
be seriously affected by the pulse nature of the incident 
wave: We find that the long-range correlation of intensity 
fluctuations increases linearly with time, 



C2(Ar,t) = C^™)(Ar)(a + /3 



-,(cw) , 
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and can exceed its cw value C2 (Ar). Remarkably, the 
above result holds for any pulse shape and duration, pro- 
vided that the correlation is measured at large enough 



times t after the extinction of the incident pulse and pro- 
vided that t is well below the Heisenberg time tn- The 
parameters a and f3 in Eq. ^ depend on the precise 
pulse shape, and (for a given pulse shape) on the ratio 
tp/to of pulse duration to the Thouless time. Besides, 
a exhibits a weak dependence on the optical thickness 
L/l of the waveguide, the dependence becoming less and 
less pronounced as L/i increases. Interestingly, we find 
that a similar linear growth with time should also be typ- 
ical for other types of mesoscopic correlations (e.g., the 
correlations of intensities corresponding to different po- 
larizations of scattered waves (23) and fluctuations (e.g., 
the fluctuations of the total transmission coefficient of a 
disordered sample jj]), provided that they are measured 
for a pulsed incident wave. Below we sketch the main 
steps of our analysis. 

Let a pulse V'o(0 be incident on the surface z = of 
a quasi- ID disordered waveguide of length L » ^ S> A 
and with cross-section A. In the diffusion approxima- 
tion, both the average intensity and the short-range cor- 
relation function of intensity fluctuations for the multi- 
ply scattered wave inside the waveguide can be found by 
summing ladder diagrams, similarly to the cw case |l3l |: 
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where (/)„(z) = (2/L)^/^sin(n7rz/L), Ar = r — ri. At = 
t — ti > 0, and we introduced 

ti 

/(")(t,ti)- /dt'e-"'(*i-*')/*°Vo(t')V^o(i' + At) (5) 



In what follows, we consider two different pulse shapes: 
a rectangular pulse |V'o(i)| = (l/ip)n(2i/ip), where 
n(a;) = 1 if |a::| < 1 and n(x) = otherwise, and a 
Gaussian pulse |'0o(t)| — {2/ y/JTtp) exp{~At'^ /tp), both 
of width tp and centered at i = 0. For these two pulse 
shapes, Eq. jSJ can be evaluated analytically, but we omit 
the exact expressions for they are rather cumbersome. In 
the present paper, we are mainly concerned with 'long 
times' that we define by requiring t 3> ip, t_D for the rect- 
angular pulse and t » t^/tDTto for the Gaussian pulse. 
At such long times, the n — \ terms give the principal 
contributions to Eqs. (jSJ and Q. We find 

if |At| < tp and /^^^(t, ti) = otherwise for the rectan- 




gular pulse, and 

ml Ifj (^^ 

for the Gaussian pulse (we introduced T = {t + ti)/2). 

As follows from Eqs. (O, © and 0, the average in- 
tensity decays exponentially with time, consistent with 
previous studies 'j, '5] . The short-range correlation func- 
tion of intensity fluctuations Q decays rapidly to zero 
for Ar > A, just as for cw illumination. Hence, the fl- 
nite duration of the incident pulse does not affect the 
short-range spatial structure of the speckle pattern. The 
correlation time tc of intensity fluctuations is set by tp 
for the Gaussian pulse, whatever its duration, whereas 
tf, ~ min(tp, t£,) for the rectangular pulse. This difference 
between long rectangular and Gaussian pulses is due to 
the intrinsic difference in their shapes: a long Gaussian 
pulse has slow variations of intensity, while the inten- 
sity of a long rectangular pulse still jumps abruptly at 
t = ±tp/2. 

Let us now consider the correlation between distant 
speckle spots, i.e. the long-range correlation of intensity 
fluctuations. This can be calculated either by joining 
four time-dependent ladder propagators using a Hikami 
box, in much the same way as it is done for continuous 
wave illumination IS] , or by generalizing the steady-state 
Langevin approach 19; 
result: 



Both methods yield the same 
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d=^r' / dt' / di'i 



V'G(r,t;r',t')-V'G(ri,ti;r',i'i) 

x(«(r',t')«(r',i'i))3hort (8) 

where G(r, i; r', t') is the Green's function of the diffusion 
equation and the spatial integral is over the volume V of 
the sample. In the rest of the paper, we restrict ourselves 
to i = ii. Substituting Eq. 10} into Eq. © we end up 
with 
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T„„p, = t-^^ / dt' / di'i ^-^n\t-t'),to-n\t-t',)/to 
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FIG. 1: Parameters a and /3 that determine the time de- 
pendence of the long-range intensity correlation function of 
intensity fluctuations C2(Ar,t) = Cj (Ar)(a -I- /St/to) in 
transmission of rectangular (dashed lines) or Gaussian (solid 
lines) wave pulses of duration tp through a quasi- ID disor- 
dered waveguide (Cj ™ is the value of correlation for contin- 
uous (cw) illumination, to is the Thouless time, and the limit 
L/l ^ CX3 is taken). The inset shows the fit to the experimen- 
tal data of Ref. |2lll (symbols) obtained by fixing L/l = 10, 
L/ La = 3 and using tp/to as a fitting parameter. 



At times exceeding tu after turning on of a very long 
{tp 3> t]j) rectangular pulse, but before its extinction, 

Eqs. ©-lini) yield the cw result C;^™^ = 2/3g. At long 
times, we can keep only the leading terms in Eqs. © 
[i.e., terms decreasing as exp{—2t/tu)] and ^ [terms 
decreasing as exp(— t/tu)]. This immediately gives Eq. 
(0) with 



/3 = 



3 smhitp/tjj) — tp/t 



D 



27r2 



smh^ {tp/2tD) 



for the rectangular pulse and 



/3 = 
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(12) 
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for the Gaussian pulse. The derivation of the above 
analytic expressions is made possible by the fact that 
only one term of Eq. ^ (the term corresponding to 
m = n=p^q = l) contributes to (3. In contrast, 
the calculation of a requires the summation of infinite 
series in Eq. 0. We have carried out this summation 
numerically for L/£ — > cx) and show the results in Fig. ^ 
together with the analytical formulas H12|) and H13|) for 
f3. 

Figure ^ suggests that two different regimes can be re- 
alized depending on the pulse duration. For short pulses 
{tp ^ i/)), we find a ex {tp/ti))^''^ and f3 ex tp/tjj for 
both pulse shapes (see Table QJ, and the long-range cor- 
relation exceeds its cw value if i > ifj/ip. In this regime. 



the incident pulse is too short for its shape to be resolved 
in transmission, and hence we anticipate that apart from 
numerical factors, our short-pulse results are fairly uni- 
versal (i.e., they hold for any pulse shape). A different 
regime is realized for pulses of long {tp ^ to) and inter- 
mediate {tp ^ to) duration: The long-range correlation 
function exceeds its cw value through the whole range of 
validity of our analysis and appears to be sensitive to the 
precise pulse shape (see Fig.^and TablcQJ. Interestingly, 
for arbitrary relation between tp and tjj the parameter j3 
seems to follow a 'universal' scaling relation (5 oc tc/to ■ 

The growth of long-range correlation of intensity fluc- 
tuations with time for Gaussian pulses of width tp '^ t^ 
has been clearly observed in experiments of Chabanov 
et al. (see Fig. 3 of Ref. |21|), thus providing an ex- 
perimental support of our result Q . Unfortunately, the 
two available data points are insufficient to deduce the 
functional dependence of C2, either on t or on tp. Never- 
theless, the data can be fit with our theory (see the inset 
of Fig. ^, provided that absorption is incorporated in 
the theoretical analysis. Absorption introduces an addi- 
tional factor exp[— (i — t')/t(j] in the expression for the in- 
tensity Green's function (here ta is the absorption time). 
As a consequence, rri^ in Eq. (|ll|l for Tmnpq is replaced 
by 771^ -|- {L/ttLo)^ (and similarly for n^, p^, and q^) 
with La = {ctaf-liY^"^ the absorption length. By fixing 
L/l = 10 and L/ La = 3 |l5Ll2ll| we obtain a satisfactory 
fit to the data with a reasonable value of tp/tu = 1.7. 

The rise of the long-range correlation function C2 with 
time can be understood physically by taking advantage 
of the path picture of wave propagation in disordered me- 
dia. In this picture, Ci is due to the crossing of two wave 
paths inside the sample 23). For a path of length s, the 
probability of crossing is given by the ratio of the vol- 
unre of a curvilinear tube of length s and diameter '-^ A 
to the total volume V = AL of the waveguide, yielding 
C2 ^ s)? /V . For cw illumination, the typical path length 
s ^ L"^ /I yields C2 ^ 1/5- For a pulsed source, s = ct 
and we obtain C2 ^ t/tu^ where tn = gto- This sim- 
ple qualitative reasoning is in complete agreement with 
the more rigorous Eq. (O since the latter can be readily 
rewritten as 



(cw) , 



C2(Ar,t) = aCr'{^r) 



rtH 



(14) 



The vanishing of a and (3 (and hence of C2) in the short 
pulse limit tp/to -^ (see Fig. ^ and Table P) can be 
understood by noting that in order to contribute to C2, 
the two crossing wave paths should arrive at the crossing 
point within a time interval ^ tp. Obviously, this be- 
comes progressively less probable as tp decreases below 
the time of wave diffusion through the disordered sample 
t]j, leading to vanishing of a, (3, and, consequently, of C2 
as tp/to tends to zero. 

Equation H14|) is expected to hold in higher- 
dimensional (2D, 3D) disordered systems as well, with 



TABLE I: Summary of asymptotic results for parameters a and /3 governing the ratio of the time-dependent long-range 



correlation function of intensity fluctuations C2(Ar, t) to its value C^ (^i") for continuous wave (cw) illumination; C-z/C^ 



t(cw) 



a + fit /to [to is the Thouless time and tp is the duration of the incident 


pulse). 
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Long pulse 


«y¥ 0.67-^1^ 
/5 ^^ ^ 




3 tp 

(27r)3/2 to 


-0-2 (^)' 

3 tp 
(27r)3/2 to 



a and (3 depending not only on tp/to but also on the di- 
mensionality and the experimental geometry. It follows 
from Eq. (|14|l that the increase of long-range correlation 
in speckle patterns with time favors bounded samples be- 
cause otherwise tn ccV -^ oo and t/tH -^ 0. 

In conclusion, we have shown that when a wave pulse 
is sent into a disordered waveguide, the speckle pattern 
formed in transmission exhibits an enhanced long-range 
spatial correlation, whereas the short-range correlation 
appears to be identical to that for continuous wave (cw) 
illumination. At long times, the long-range correlation 
grows linearly with time, whatever the shape and the du- 
ration of the incident pulse, and can exceed its cw value. 
Our analysis breaks down beyond the Heisenberg time, 
where the speckle pattern is dominated by anomalously 
localized states and the diffusion picture of wave prop- 
agation is not valid any more. It is worthwhile to note 
that other types of mesoscopic phenomena should exhibit 
a similar linear enhancement with time, when observed 
for a pulsed source of waves. In particular, correlations 
between waves of different polarizations and full prob- 
ability distribution functions of transmitted intensities 
have been recently measured in dynamic microwave ex- 
periments |21| . Our analysis suggests that the theoretical 
description of these measurements at times much shorter 
than the Heisenberg time can be obtained by simply re- 
placing the dimensionless conductance g by g/{a+(3t/t£i) 
in formulas corresponding to the stationary (cw) regime. 
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